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We address the problem of stitching together the vacuum, static, planar-symmetric Taub 
spacetime and the flat Friedmann-Robertson-Walker cosmology using the Israel thin-shell 
formalism. The joining of Taub and FRW spacetimes is reminiscent of the Oppenheimer- 
Snyder collapse used in modeling the formation of a singularity from a collapsing spherical 
ball of dust. A possible mechanism for the formation of a planar singularity is provided. 
It is hoped that tackling such example will improve our intuition on planar-symmetric 
systems in Finstein’s general relativity. 
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1. Introduction 

Similar concepts that appear in different fields in physics are not uncommon. We 
know from classical electromagnetism (EM) that non-overlapping regions of space 
must be described by piecewise continuous scalar and vector potentials, but the 
derivatives of these potentials, i.e. electric and magnetic fields, need not be con¬ 
tinuous at the junctions. Discontinuities in the electric and magnetic fields at the 
junctions give rise to surface charge and current distributions. The same scenario 
can be seen in classical general relativity (GR) when one tries to patch together non¬ 
overlapping regions of spacetime described by different geometries. In this case, the 
components of the metric tensor are analogous to the 4-potential in classical EM 
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and discontinuities in the derivatives of the metric tensor give rise to surface energy- 
momentum distributions at the junctions of the spacetimes stitched together. 

A suitable formalism that can be used to patch together spacetimes of different 
geometries has been formulated by IsraeP^ for the case of non-null hypersurface 
junctions, and by Barrabes and Israel^ for the case of null hypersurface junctions. 
For the case of non-null hypersurface junctions, the conditions for a continuous and 
smooth transition at the junction between two spacetimes are: 

[hab] = 0 ( 1 ) 

[Kab] = 0 ( 2 ) 

where hab and Kab are the induced metric and extrinsic curvature on the hypersur¬ 
face, respectively. In the event that ([2|) is not satisfied, a surface energy-momentum 
distribution exists at the surface junction: 

Sab = —^i\Kab]-[K]hab) (3) 

OTT 

where K is the trace of the extrinsic curvature, and e takes the value —1 for a 
timelike hypersurface and -|-1 for a spacelike hypersurface. Meanwhile, we also have 
a similar set of junction conditions for null hypersurface junctions: 

V^ab] = 0 (4) 

[Cab] = 0 (5) 

where pab is the induced metric of the two-dimensional hypersurface and Cab is the 
transverse curvature, which is the analog of the extrinsic curvature for null hyper¬ 
surfaces. If there is a jump in the transverse curvature across the null hypersurface, 
then there exist surface energy, current, and pressure at the junction. 

Israel first used the junction conditions in joining together Schwazrschild and 
Minkowski spacetimes using a spherical thin shell of pressureless dust. It was found 
out that: (I) the rest mass of the thin shell is constant throught its history, and 
(2) the total energy of the shell may be broken down into its kinetic energy and 
binding energy.® Meanwhile, Barrabes and Israel extended the previous case to a 
null shell joining an exterior Schwazrschild region and an interior region described 
by a spherically symmetric spacetime containing a false vacuum. They showed that 
the expansion of the lightlike shell requires that the energy density of the shell is 
consumed in the process.® 

In the examples mentioned previously, interesting physical interpretations have 
been drawn in studying those systems. Now, we wish to follow a similar track in 
studying thin shells in the context of planar symmetric spacetimes. Interest in the 
physical application of planar symmetric systems is rooted from modeling domain 
walls, which are topological defects that have been formed in spontaneous sym¬ 
metry breaking processes during cosmological phase transitions. In the thin wall 
approximation, domain walls may be described by an energy-momentum tensor 
proportional to a Dirac delta.® 
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In this paper, we attempt to consistently patch together the Taub planar sym- 
metric spacetime with a flat Friedmann-Robertson-Walker (FRW) model. This 
study is motivated by extending the work of Oppenheimer and Snyder in studying 
the eventual collapse of a spherical ball of dust to a Schwarzschild singularity; we 
wish to provide a similar mechanism that might explain the formation of planar 
singularities. 

This paper is organized as follows. Chapter 2 provides a brief account of planar 
symmetric spacetimes in GR. Chapter 3 contains the main results of this paper; 
we have demonstrated how to stitch together Taub and FRW spacetimes using 
subluminal dynamic thin shells. We also looked at certain configurations of dynamic 
thin shells that will eventually lead to an eternal expansion or collapse of the FRW 
spacetime. The last chapter contains the conclusions. 


2. The Taub planar symmetric spacetime 


It was shown by TaulP that the metric for a vacuum planar-symmetric spacetime 
with zero cosmological constant can be written as: 

= (—-I- dz^) + (1 + kz) {dx^ -\- dy'^) (6) 

v \ k z 


By suitable coordinate transformations, the above may be recast into the following 
formsPI^ 


ds^ 


-f = [-dt^ + d(^) + ( [dx^ + dy'^) , C > 0 
vC 


( 7 ) 


ds^ =- 0 ( 8 ) 

One can immediately notice that the Taub spacetime is static since the metric is 
invariant under translations along the t coordinate. Also, by looking at ([8]), one 
may ascertain that there is a singularity at z = 0. To show that this is a physical 
singularity, we may evaluate the Kretschmann scalar, and we And that: 

K ^ ^ (9) 

This curvature invariant clearly blows up at z = 0, and thus signifies that z = 0 
is indeed a physical singularity. The hypersurface z = 0 has been interpreted by 
Bedran et. al. as an infinite mass sheet with negative energy density and is repulsive 
in nature.^ One may perform a quick calculation to verify the repulsive nature of 
this planar singularity. We consider an observer that wishes to be fixed at z = zq. 
For simplicity, we neglect the motion along the x and y directions. In this case, we 
may write down the four-velocity of such an observer: 

= \zo\'^^^dt 


The acceleration of this observer in the coordinate basis is: 

1 


( 10 ) 
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Looking at one can see that for an observer to stay fixed at z = zq > 0, the 
observer must maintain an acceleration pointing towards the planar singularity at 
z = 0, indicating that the planar singularity is repulsive. Similarly, for z = zq < 0, 
the observer must maintain an acceleration pointing towards the planar singularity, 
and we again arrive at the conclusion that the planar singularity is indeed repulsive. 

3. Stitching Taub and FRW spacetimes 

At this point, we want to stitch together Taub and Friedmann-Robertson-Walker 
(FRW) spacetimes, motivated by the extension of the work of Oppenheimer and 
Snyder in studying a spherically symmetric collapse of pressureless dustP to planar 
symmetric collapse. In this spherical collapse, the exterior geometry is Schwazrschild 
while the interior geometry is that of a closed FRW spacetime. A modern treatment 
by PoissorP showed that a smooth boundary between the two spacetimes implies 
that the mass of the star M can be related to the energy density p of the FRW fluid 
via: 

M = (12) 

O 

where R is the radius of the ball of dust. 

We have chosen to work with a fiat FRW spacetime, whose metric is given by: 

ds^ = -dtl_ + a^(t_) (dC^ + dx^ + dy^) (13) 

where a{t-) is known as the scale factor, which accounts for the homogeneous 
expansion (or collapse) of a Euclidean 3-volume. 

We shall use the Taub metric in describing a static planar symmetric vacuum 
spacetime: 

ds^ = -+ dz"^ + {dx^ + dy^) (14) 

3.1. Description of the Planar Junction 

Consider the hypersurface junction E dividing spacetime into two regions: region I 
(FRW) and region II (Taub). Region I covers 0 < ^ < Z_(t) and region II covers 


z > Z+(t). As seen from region II, E can be parametrized by t+ = T+(t) and 
z = Z+(t). The tangent basis vectors on E as seen from region II are: 

= f+dt^ + Z+d, (15) 

epd^ = d, (16) 

= dy (17) 

The overdot stands for differentiation with respect to t. Hence, the induced metric 
on E as seen from region II is: 

ds| = dr'^ + Z+/^(dx2 -h dy"^) (18) 
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Meanwhile, as seen from the FRW side, the junction is parametrized by t- = T-{t) 
and C = Thus, the tangent basis vectors are: 

u-°‘d^ = f_dt_ + Z_d^ (19) 

e~°'da = ( 20 ) 

= dy ( 21 ) 

so that the induced metric on S as seen from region I becomes: 

dsl = -(t^- a^{T_)Z'Dj + a'^{T_){dx'^ + dy"^) (22) 

The first junction condition and the normalization condition for 4-velocities 
imply the following relations: 

T! - a^{T_)Z^ = ZZf^'^ (tI - zf^Zl^ = 1 (23) 

zj^ = a2(r_) (24) 


4. Calculations Regarding the Second Junction Condition 


Let us now calculate the extrinsic curvature on each planar junction. Before we 
begin, we note that the normal one-form on S as seen from the FRW side is: 


n^dx°' = 

|a(T_)| -f T_dc) 

(25) 

and from the Taub side: 

|Z+|i/3 ( Z+dt++T+dz) 


P 

II 

(26) 


Note that we have chosen the normal vector to point from the FRW side to the 
Taub side, i.e. > 0 and n‘’ > 0. The extrinsic curvature can be computed via: 

Kab = (27) 


Let us begin calculating the extrinsic curvature on the Taub side. We have: 


Kt-r = n'.oU ' “u 


a\p 


= -n+a+“ 


(28) 

(29) 


where is the 4-acceleration of an observer comoving with the thin shell in the 
Taub side. However, 


a+‘ = 


dr 

= T+- 


a+^ = 


du+^ 
dr 

= Z+- 


+ r+fu+‘^u+d 

(30) 

^ T+Z+ 

3Z+ + + 

(31) 


(32) 

1 rp2 

3Z+/^ 

(33) 


and: 
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Thus, using (1^ . (IHTl) . and (1551) on (1551) . we have: 


3z:; 

wo 


|Z+|l/3 


Z+T+ - T^Z. 


+ ^+ 


(34) 


Meanwhile, the other non-vanishing components of the extrinsic curvature are: 

= K+y = ^sgn(Z+)r+ (35) 

We now move on to the FRW side. Similar to the previous calculation of 
in the Taub side, we first compute the components of the 4-acceleration of the thin 
shell in the FRW side. We have: 


a"* = T_ 


a-^ = Z_ 


a{T-) da{T) 


T_ dr 
2 da{T- 


Zt 

Z_ 


(36) 

(37) 


a(T_) dr 

Thus, we have the following component of the extrinsic curvature in the FRW side: 

K-^ = |a(r_)| (z.r. - 


-sgn(a(r_)) 


da(T_) ( Z_ 


dr 


r_ 


(i+ p) 


Also, we have: 


Kxx = K^y = W(T_)sgn(a(T_)) 


(38) 


(39) 


4.1. Is a boundary junction possible? 


In the Oppenheimer-Snyder collapse, the surface of the spherical fluid of dust is 
a boundary layer; a smooth junction exists between the FRW interior and the 
Schwarzschild exterior. We also ask if this is also the case in stitching Taub and 
FRW spacetimes. 

For simplicity, we assnme that the junction traces a geodesic in either spacetime. 
Hence, in the Taub side, the 4-velocity is: 

= E ^ (40) 

where E is a constant of motion due to a translation invariance of the Taub 
metric. In the FRW spacetime, we have: 


Uql — 


v\ 


\TJ) 


+1 a*. + 


Pc 


\TJ) 


(41) 
















September 29, 2015 0:43 WSPC/INSTRUCTION FILE ws-ijmpd 


Dynamics of a planar thin shell at a Taub-FRW junction 7 


where pq is a constant of motion due to a C translation invariance of the FRW metric. 
Since the trajectory of the junction on either spacetime is a geodesic, comoving 
observers with the shell experience zero acceleration, and thus iLjlj. = K~^ = 0. 
Now, plugging in the components of the 4-velocity of the shell to (l3^ and (l3^ and 
invoking the second junction condition, we arrive at: 


PC sgn(a(r_)) da(T_) 



(42) 



,2 I n2(T \ dr 


Using (l24l) . we have a differential equation for a(r_): 


PC da(T) 



(43) 



dr 


One can show, by using (H51) and dnj, that the FRW scale factor turns out to be: 



(44) 


We have arrived at the conclusion that a smooth junction, where the junction 
follows geodesic motion in either Taub or FRW spacetime, must have a scale factor 
given by (1441) at the FRW side. It might be interesting to probe if the system 
of equations, obtained by invoking the first and second junction conditions, has 
a unique solution up to certain arbitrary constants, which may be obtained upon 
setting initial conditions. 

4.2. Friedmann equations in the presence of a boundary layer 

The Friedmann equations may be obtained by plugging in the FRW metric into the 
Einstein field equations and assuming a cosmological perfect fluid that permeates 
in FRW spacetime. The Friedmann equations for a flat FRW model are: 



(45) 


(46) 


where p is the density of the cosmological fluid, P is the pressure of the cosmological 
fluid, and the overdots in (1451) and (l46l) stand for differentiation with respect to . 
Plugging in (|44[) into the Friedmann equations will yield: 




(47) 



( 48 ) 


2E 


where: 
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The equation of state p = —P for this FRW fluid corresponds to a vacuum- 
dominated FRW universeP 

In cosmology, when one models the nniverse as a flat FRW spacetime contain¬ 
ing radiation, matter, and vacuum energy, the evolution of the universe can be 
marked by different epochs in its expansion: radiation-dominated at early times, 
then matter-dominated, and finally a vacuum-dominated expansion!^ During a 
vacuum-dominated expansion, the scale factor is exponentially increasing with FRW 
coordinate time, i.e. a{t) oc and hence ingoing/out going photons emitted at 
some coordinate time ta > 0 will only reach a spatial distance of: 


d-H oc 


dt' 

a{t') 


= H^ 


^-Hta 


(50) 


Photons cannot reach beyond this spatial distance. This might lead us to suspect 
that somewhere beyond the cosmological horizon, there exists a boundary layer that 
separates the FRW spacetime from the Taub spacetime, but cannot be observed 
from our current location. 


4.3. Stationary thin shell in FRW 

Let us study the case where the junction is stationary in FRW spacetime, i.e. Z_ = 
0. Even by considering this simple example, it will turn out that we can extract 
interesting results such as qualitative descriptions regarding the trajectory of the 
thin shell. Looking at (IMl) and (IMl) . we see that the components of the extrinsic 
curvature of the hypersurface junction as seen from the FRW side vanish, i.e.: 


Kb = 0 (51) 

Meanwhile, the components of the extrinsic curvature in the Taub side, ([34)l and 
((35)l . will stay as they are. Notice that unlike the previous case of the thin shell 
tracing a geodesic trajectory, it appears that the Taub and FRW spacetimes cannot 
be joined smoothly in this case because if we set [Kab] = 0, we will have T+ = 0, 
which will make the 4-velocity of the planar shell to be spacelike in the Taub side. 
Thus, the planar sheet must be made up of matter whose surface energy-momentum 
tensor described by: 


= - — {[Kab] - [K]hab) 

OTT 

= -^{Kb-Khab) 


(52) 

(53) 


Sab = - 


67rz; 


:473 


T+ 

247r 

Stt 


0 


(zl 


p 1 

Z+T+ - T+Z+ 

0 


T+ 

247r 


Sr (2 


{zl + 1 ) 


Z+T+-T+Z. 


+ ^+ 


(54) 
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If the planar shell is made up of a perfect fluid, then: 

Sab = 


cr 0 0 

0 pa^ 0 
_ 0 0 pa^_ 


(55) 


where a is the surface energy density of the shell and p is the pressure on the shell. 
Hence, we have: 




6ttZ_ 


4/3 


(56) 




-p = -- 


247rZ: 


4/3 


{zl + i) 


8ttZ. 


1/3 


{z+f. 


T+ - T+Z+ 


(57) 


We may rewrite our expressions for the surface energy density and surface pressure 
in terms of the scale factor a(r). From (1231) and (IMl) . and noting that the thin shell 
is restricted to be stationary with respect to the interior FRW region, we have: 


Z+ = -a^/^{T)a{T) 

t+ = a}/‘^{T) + ^a(T)d2(r)^ 


1/2 


(58) 

(59) 


and hence we eventually obtain the surface energy and pressure in terms of the 
FRW scale factor: 


cr = — 


(1 


9 ■ 2 \ 


-p = 


(l + |ad^ + |a^d) 
247ra3/2 + 


(60) 

(61) 


4.4. Consistency Check 

Now that we have the expressions for the surface energy density and pressure, 
we wish to check if our expressions for the surface energy density and pressure 
are correct. The surface energy density and pressure must satisfy the continuity 
equation: 

S% = -e[r] (62) 

where ja = Tf^pe^n^ is a surface current, and the bar denotes covariant differentia¬ 
tion on the surface, using a connection that is compatible with the induced metric 
on the junction. The reason why a and p must satisfy is because imposing the 
junction conditions on the junction between the Taub and FRW spacetimes (which 
are solutions of the EFE) implies that the EFE are satisfied everywhere and, as a 
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result, the energy momentum tensor satisfies the integrability condition = 0 
everywhere. 

Owing to the fact that the Taub spacetime is a vacuum solution and that the 
energy-momentum tensor of the FRW spacetime is diagonal, we have [j“] = 0. 
Carrying out the calculations, we have: 


oab _ 

J IL — 


1 


= 0 


a2(r)' 

Taking the a = t component of the equation, we have: 


1 d 


^(r) dr 


Eventually, we obtain: 


[a^{T)s^^] +Ti,s^^ + riySyy = o 


^ (aa^) +p^(a^)=0 


dr 


(63) 


(64) 


(65) 


Notice that (1551) has a resemblance with the first law of thermodynamics. To push 
this resemblance further, we note that the total area of the planar shell is: 


J J coiaclhrrl^^^dxdy 

(66) 

a?{T)L‘^ 

(67) 


Again, we have introduced a cut-off in place of a divergent integral over all x 
and y. Thus, we identify the “total energy” of the shell to he U = aa?LF' and the 
“total area” of the shell to be A = 

It turns out that the expressions for the energy density and pressure, given by 
dsni) and (l6T]) . respectively, satisfy the conservation relation (1651) which is to be 
expected. 


4.5. Equation of state: cr — —wp 

Let us consider an equation of state of the form: 


cr = —wp 

for the dynamic thin shell, where w 0. Using dsni) and (El), we have: 


(I + fad^) 


9„a2\1/2 


w 


0/0 — 1 /o I 1 H” —Qj^d 

6^3/2 247ra3/2(l + |ad2)^^" V 2 2 


9 


9 


Simplifying yields: 


9 9 

(4 — w) -I- -(2 — w)ad^ — -wa^'d = 0 


( 68 ) 


(69) 


(70) 


We may rewrite the above differential equation in terms of which amounts to 
the equation of motion of the shell in the spatial direction. Erom (l24l) , again noting 
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that the thin shell junction is stationary relative to the interior FRW region, we 
now have: 


(4 -w)(l + - 3wZ+Z+ = 0 (71) 

A straightforward integration of the above expression yields: 

2 ( 4 -») 

3 zu 

= -1 (72) 



At this point, we may now give qualitative descriptions to the fate of a stationary 
thin shell junction in the FRW region for any given w. One can view the above 
relation as an effective energy equation with an effective potential: 


V(Z+) 



(73) 


where tq is an integration constant, and energy equal to —1. For any w ^ 0, there 
is always a turning point, and hence one can identify the classically allowed and 
forbidden regions. If the exponent of the effective potential is positive, i.e. 0 < w; < 4, 
then the thin shell avoids an eventual collapse since the classically allowed region is 
[ro,oo), where tq is the location of the turning point. On the other hand, if ic > 4 
or w < 0, then the classically allowed region is (0,to], and hence the thin shell will 
eventually collapse. 


4.6. Special Case: Domain wall junction 


Supposing that the planar junction is a domain wall, we impose the equation of 
state cr = — ppdlllll A domain wall would correspond to the case w = 1. From 
the previous section, we then see that this domain wall avoids an eventual collapse, 
and thus the FRW region will expand forever. From (I72|) . we may write down the 
equation of motion of the domain wall: 


Setting ^+(0) = Zq, we now have a solution for Z+{t): 


Z+(t) = Z+o cosh- 1-770 


(74) 


(75) 


where 770 = cosh ^ • We can also rewrite the solution as: 

Z+{t) = Zq cosh -I- ^Jz^-T^ sinh (76) 

Looking at our expression for Z^{t) reinforces our qualitative analysis on the dy¬ 
namics of the thin shell in the previous section: the thin shell avoids an eventual 
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collapse to a singularity. Meanwhile, the expression for the scale factor a(T) will 
immediately follow: 


i(r) = 


cosh + \JZq — Tq sinh ( — 


2/3 


(77) 


V '^0 / '' V t-q , 

Finally, the surface energy density and pressure of the planar junction are given by: 


1 


a = —p = 


Gttto 


1 - - tanh^ (— + Vo 
9 \to 


1 1/2 


(78) 


4.7. Friedmann equations in the presence of a domain wall 
junction 

Writing: 


a(T) = T^ 


^/^osh2/3 


To 


+ m 


and plugging into the Friedmann equations will yield: 

1 


P = 




tanh 


To 


Vo 


P = - 


1 


2TrT§ 


(79) 

(80) 
(81) 


Rewriting the energy density of the cosmological fluid so that it contains the pressure 
of the fluid: 


p = -iptanh^ ( — + Vo 
3 V To 


(82) 


To aid us in understanding the physical situation, we note that the z-coordinate of 
the planar junction decreases from infinity to z = Zq in the proper time interval 
(—oo, —ryoTo), and then bounces off from z = tq to infinity in the proper time 
interval (—ryoTo,c»). In this case, we have an initially contracting FRW spacetime, 
and then expands after the planar junction reaches the bounce location. Meanwhile, 
we can see that the density of cosmological fluid in the FRW side decreases from 
to zero when the planar junction moves closer to the bounce location, and 
then increases again asymptotically to when the planar junction has bounced 
off z = Tq. 

It can be shown that in the presence of this domain wall junction, the FRW 
spacetime possesses a cosmological horizon during its expansion phase. Begin by 
considering the integral 


dn = 


dr 

a(r) 


(83) 
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We require dn to be finite so that this FRW spacetime possesses a cosmological 
horizon. Now, 


dn = Tr 


-2/d. 


dr 


Consider the integral: 


_t1/3 
— ^0 


h cosh^/^ +770) 
du 


sh^/^ 


< T, 


1/3 


Irio cosn ' u 
du 


1-00 cos 


h2/3 


u 


- .^1/3 


I = T, 


du 


1-00 cos. 




(84) 

(85) 

( 86 ) 


(87) 


Once we have shown I to be finite, it then follows that dn is bounded above, and 
hence I tends to a finite value. Letting x = e“ gives: 

^ ^ . 2 ^ 112 / 3 ^^ 

Jo (X^ + 1) ' 

Now, let X = tan0. Thus: 


/ = 



(89) 

(90) 


In the last line, we used the integral definition of the Beta function: 

t/2 


B{m + I, n + I) = 2 / 
Jo 

Hence, we have the inequality: 

0 < d,, < (^) 


Cos2m+lg,sin2n+l0^0 




(91) 


(92) 


and thus we have shown that there is a cosmological horizon in this FRW spacetime. 
Hence, there will come a time when the domain wall will be obscured from view 
with respect to observers sufficiently far away from the wall. 


5. Conclusions 

We have been able to demonstrate that one can use the Israel thin shell formalism 
to join together the Taub and FRW spacetimes, and deduce some physical conse¬ 
quences. In stitching the Taub and flat FRW spacetimes, we first obtained expres¬ 
sions for the extrinsic curvature of a generic spacelike hypersurface junction joining 
these two spacetimes. We have shown that, similar to the Oppenheimer-Snyder 
collapse, one can have a smooth transition at the junction of these spacetimes. 
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Meanwhile, for a junction that is stationary with respect to the FRW region, the 
Taub and FRW spacetimes cannot be joined smoothly. As a consequence of the sec¬ 
ond junction condition, the junction between these spacetimes possesses a surface 
energy density a and surface tension —p, which are dependent on the proper time 
of an observer comoving with the planar junction. We then considered an equation 
of state a = —wp, and gave qualitative descriptions on the motion of the thin shell 
for two distinct cases: 0 < ic < 4 (case 1) and w < 0, ic > 4 (case 2). Case 1 
does not form a planar singularity, while case 2 eventually leads to a collapse of 
the thin shell junction to a singularity. We also considered the special case w = 1, 
which corresponds to a domain wall, and obtained an analytic expression for the 
scale factor and the spatial location of the domain wall. It has been shown that 
the domain wall cannot form a singularity, which is consistent with our qualitative 
investigations. 

A natural extension to this work is to look for a possible relationship between 
the equation of state of the surface energy and pressure of the Taub-FRW thin 
shell junction and the the equation of state of the cosmological fluid. One could 
also consider general forms for the scale factor a^r) and deduce the properties of 
the thin shell junction. Finally, one could also consider stitching planar symmetric 
spacetimes other than the vacuum solutions to the field equations. 
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